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1.  INTRODUCTION 

Research  results  In  various  areas  such  as  speech  processing  [1], 
array  processing  [2],  and  adaptive  tracking  [3]  have  indicated  that  the 
gradient-lattice  algorithms  have  superior  convergence  behavior  over 
classic  time-recursive  methods  such  as  Widrov's  LMS  approach  [4].  The 
majority  of  this  work,  however,  has  been  of  a  primarily  experimental 
nature.  It  is  the  purpose  of  this  paper  to  provide  a  theoretical  frame¬ 
work  for  the  convergence  study  of  such  gradient-lattice  algorithms.  Since 
the  convergence  of  adaptive  algorithms  is  primarily  a  stochastic  problem, 
a  mean-square  criterion  will  be  considered  as  opposed  to  the  convergence 
of  the  mean  of  the  reflection  coefficients,  which  is  mostly  used  in  the 
literature  [1,3,5-10]. 

A  stochastic  version  of  a  fixed-point  theorem,  which  was  developed  by  • 

Oza  [11-12],  will  provide  the  basis  in  establishing  convergence  conditions 
for  the  gradient-lattice  algorithm.  Such  an  approach  has  been  taken  in  an 
earlier  paper  [13]  to  study  convergence  properties  of  the  LMS  algorithm 
and  is  based  on  the  notion  of  a  contraction  mapping  on  a  "stochastic” 

.1 

Hilbert  space  [13].  In  the  same  framework,  expressions  for  the  convergence  ! 

J 

rate,  conditions  on  the  stepsize,  and  misadjustment  are  derived  by  con- 

j 

J 

sidering  a  "distance"  measure  in  the  appropriate  space.  j 

II.  MATHEMATICAL  BACKGROUND 
1.  The  Lattice  Filter 

The  lattice  filter,  as  depicted  in  Figure  1,  will  formally  be  con- 

* 

sidered  as  the  cascade  of  elementary  operators  A^  (Fig.  2): 

N 

A  -  ir  A  (1) 

i-1 


t 
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where  N  is  the  order  of  the  filter.  Since  the  input  of  each  element  in 
the  lattice  can  be  represented  by  the  "augmented"  vector: 


fi_l(n) 

-bi-l(n)- 


(2) 


A^  will  be  assumed  to  be  an  operator  on  >  the  Hilbert  space  of  zero- 

mean,  wide-sense  stationary  stochastic  (2x1)  sequences  with  inner  product: 

<8^0 »  ej(*)>  *  E{e^(n)  e^n)}  (3) 

It  is  also  assumed  that  e^(<)  has  finite  average  power*.  The  correlation 
matrix  at  the  input  of  each  stage  is  given  by: 

E{fii(n)  §i(n)}  " 


where,  by  symmetry,  the  average  power  of  f^n)  and  b^(n)  (E^)  are  equal, 
can  be  written  in  matrix  form  as: 


S^n) 


1  k^z 


-1 


-1 


l_ki  z 


1  0 

0  z"1 


Si_l<n> 


1  k. 


ki  1 


s  u  ^  a1_1(n) 


(5) 


it 

The  general  theory  developed  in  this  paper  is  applicable  to  non- 
Gaussian  stochastic  processes.  However,  simple  expressions  for  the 
appropriate  norms  can  be  obtained  if  Gaussian  assumption  is  made. 


t 
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where  z”1  is  the  unit  shift  operator,  and  U  and  are  the  matrices 
defined  by: 


where  the  magnitudes  of  the  "reflection  coefficients"  are  always  taken 
to  be  smaller  than  one.  The  norm  of  [14]  is  given  by: 

llAjI  -  {|lUKifii_1(*>l|}  <7> 

llfii_l(*)  11  Ml 

where  ||  •  ||  denotes  the  norm  induced  by  the  inner  product  in  (Eq.  3). 

Since  it  is  straightforward  to  show  that  U  is  a  unitary  operator,  (Eq.  7) 
reduces  to  the  simple  matrix  norm  of  K^,  which  is  given  by  its  spectral 
radius  [14]: 

II  A,  II  -  I’Wxl  '  1+lkil  <8*> 

where  ^max  is  the  largest  eigenvalue  of  A  lower  bound  on  the  output 

power  can  be  obtained  by  replacing  sup  by  inf  and  using  the  same  proce¬ 
dure.  This  gives: 

U-|ki|)||ei_1(*)||  <  ll§i(->||  <  (i+|ki|)||ei_1(-)||  (8b) 

Note  [10]  that  the  lower  bound  is  actually  attained  when  e^_^(»)  consists 
of  the  actual  innovations  processes  (optimal  linear  predictor).  (Eq.  8b) 
was  previously  obtained  by  Makhoul  using  different  methods  [10]. 


i 
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2.  The  Adaptive  Lattice  Filter 

Thus  far,  we  have  considered  the  lattice  filter  as  a  fixed,  deter¬ 
ministic  operator  on  the  space  H2(fi).  If  the  lattice  is  used  in  an 
adaptive  algorithm,  however,  two  major  differences  arise.  First  of  all, 
if  a  time- recursive  method  is  used,  the  operators  A^  will  become  time- 
varying  and  thus  the  assumption  of  wide-sense  stationarity  of  the  input 
is  no  longer  valid.  In  the  sequel,  however,  we  will  assume  that  the 
stepsize  in  the  adaptive  algorithm  is  sufficiently  small  such  that  the 
inputs  to  successive  stages  are  at  least  locally  stationary.  Secondly, 
the  parameters  determining  the  operators  A^  will  be  based  on  measured 
data,  thus  resulting  in  a  stochastic  operator.  Consequently,  the 
assumption  of  Gaussian  statistics  of  the  input  to  each  stage  will  be 
violated.  Again,  under  the  same  assumption  of  a  small  stepsize,  it  can 
be  assumed  that  the  statistics  of  the  sequences  e^C*)  will  be  suffi¬ 
ciently  Gaussian. 

In  order  to  calculate  the  mean-square  norm  of  the  stochastic  lattice 
element,  we  can  again  use  (Eq.  5)  where  now  the  reflection  coefficients 
in  the  matrix  are  random  variables: 


'Si-1 


sup 

<•>11-1 


<I|K 


e  ( • ) 
i  si-lv  ’ 


sup 
(*) 


E{e^_i (n)  e^^n)} 


‘i-1 


Under  the  widely  used  assumption  that  the  operator  and  the  sequence 
^(*)  are  uncorrelated,  (Eq.  9)  becomes: 


i 


(9) 
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|| A  || 2  =  sup  E{eJ  (n)  E{K2}  e  (n)} 

1  He1-1(')ll  ”1 

<  ||e{K2})|  (10) 


(by  the  Schwarz  inequality  and  the  definition  of  matrix  norm). 

Because  the  matrix  norm  is  attained  by  an  eigenvector  and  since  in  this 
case  the  equality  sign  in  the  Schwarz  inequality  holds,  we  have: 

HVi2  -  iie(k2)ii  do 


Noting  that: 


EO^}  =  E  { 


’"l+k2 
1 


2k. 


l+cr2+k2  2k . 

ii  l 

2-2 

2k  1+0  +k, 

i  l  i 


where  a ^  and  are  the  variance  and  mean,  respectively,  of  k^,  we  can 
write  bounds  for  1 1 e ^ ( • ) 1 1  similarly  to  (Eq.  8b): 


(12) 


[a2+(l-|£i|)2]||ei_1(.)||2  <  He^-)!!2  <  [o2+(l+|C.|)2]||ei_1(-)||2  (13) 


3.  The  Stochastic  Fixed-Point  Theorem 

The  following  theorem,  which  forms  the  basis  for  the  discussions  in 
the  sequel,  is  a  stochastic  fixed-point  theorem  and  was  used  by  Oza  [11,12] 
in  a  system  identification  problem.  For  the  proof,  we  refer  to  [11]. 
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Theorem 

.  GO 

Let  {T  a  se9uence  of  (random)  operators  on  a  Hilbert  space 

H  and  let  Tn-*-T  where  T  is  a  contraction  mapping,  i.e. 

iim  ||Tny-T  y||  =  0  V  y  e  H 

n-wo 


and 


l|Ty1-Ty2||  <  Ily1-y2||  V  y1  y2  e  H 

Then  the  sequence  generated  by: 

yn+l  =  Tn  yn 

yg  fixed,  but  arbitrary  in  H,  converges  strongly  to  the  fixed-point  of  T. 

III.  CONVERGENCE  PROPERTIES 

1.  The  Adaptive  Time-Recursive  Algorithm 

It  is  well-known  [15]  that  the  lattice  filter  is  a  natural  implemen¬ 
tation  of  the  Levinson  algorithm  to  solve  the  linear  prediction  problem. 
In  the  "optimal"  case,  the  reflection  coefficients  k^  are  determined  in 
the  Levinson  recursion  by: 

Ei-1  ki  *  _  Ci-1 

and  the  sequences  f^(*)  and  b^(*)  are  the  innovations  processes  of  the 
ith-order  forward  and  backward  prediction  [16],  Practically,  however, 
neither  the  average  power  (E^)  nor  the  crosscorrelation  (C^)  of  the 
innovations  processes  are  completely  known,  and  an  adaptive  method  then 
consists  of  calculating  the  reflection  coefficients  based  on  estimates 
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of  Cj,  and  which  are  obtained  from  actual  data.  A  very  common  and 
simple  form  of  adaptation  is  achieved  when  the  average  values  in  (Eq.  4) 
are  replaced  by  their  instantaneous  values  and  the  solutions  to  (Eq.  16) 
are  obtained  time-recursively.  This  leads  to  an  adaptive  algorithm  given 
by  (fr  Ref.  [3]): 


i,n+l 


i,n 


-  a 


IEl-l<n)  ki,n  + 


(17) 


where 


ei-l(n>  “  I  [f^OO+bj^OO] 


and 


(18) 


£i-l(n)  '  fi-l(n)  bi-l(n)  (19) 

Note  that  if  the  estimates  are  taken  to  be  the  exact  values  (i.e., 

£j_^(n)  ~  C^_^(n)  * Cj_^) ,  the  sequence  in  (Eq.  17)  converges  if 

2 

<  g  •  the  same  condition  as  is  obtained  for  convergence  of  the  mean 
of  the  reflection  coefficients  [3]. 

2.  Conditions  on  the  Stepsize  (a) 

Strictly  speaking,  an  adaptive  lattice  element  A.  converges  to  a 

1 9n 

lattice  element  (the  "optimal"  operator)  if: 

lim  ||Ai  n  ii_1(*)  8^(0  ||  "  0  (20) 

n-Kx>  ’ 

Practically,  however,  most  adaptive  algorithms  can  only  result  in  out¬ 
puts  which  are  within  a  certain  "distance"  of  the  optimal  value,  i.e.. 


lim 

n-*» 


II  A, 

11  i,n 


si-l^  ~  Ai 


®i-l ( *  > 


M. 


<  OD 


(21) 


l 
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where  is  the  (unnormalized)  misadjustment .  Using  (Eq.  5),  the  norm 
in  (Eq.  21)  can  be  written  as: 


-Ai>  Si-jO 


(22) 


E{lki,n-kil>  E 


i-1 


(23) 


Thus,  convergence  of  the  adaptive  element  is  determined  by  convergence 
of  the  reflection  coefficients. 

In  order  to  obtain  limits  for  the  stepsize  a,  (Eq.  17)  can  be  written 
in  terms  of  a  stochastic  fixed-point  theorem  as  follows:  write  (Eq.  17)  as: 


ki,n+l  “  ki  ‘  (ki,n  "  V  [1  ~  “  ^i-l(n)  1 

-  «tCi_i(n)  +ki  §i  l(n)] 


(24) 


and  define  the  stochastic  operator 
variables  by: 


T.  on  the  space  of  Gaussian  random 
i.n 


Ti>n  k  =  [l-aE1_1(n)]k-  a[Ci_1(n)+ki  E^_^(n) ] 
for  all  k.  Then,  (Eq.  24)  can  be  expressed  as  a  fixed-point  problem: 


(25) 


(ki,n+l_ki)  “Ti,n  ^i.n-11!5 


(26) 


and  using  the  contraction  mapping  principle  of  Section  1.3,  the  recursion 
in  (Eq.  26)  will  converge  if: 

||(1-aEi_1(n)](y1-y2)||  <  ||y1-y2||  (27) 


for  y^  and  y2  fixed  but  arbitrary. 
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Since  the  norm  In  (Eq.  27)  is  given  by: 

llU -a  Ei_1(n)](y1-y2)||2 

-  E{|l-a  Ei  ;i(n)|2}  ||y1-y2||2  (28) 

this  condition  becomes  (see  Appendix  A) : 

l  +  a2[2  E2_1  +  C2_1]-2aE1_1<l  (29) 


or 


2E 


0  <  a  < 


i-1 


2  2 
2Ei-l  +  Ci-l 


(30) 


Various  interpretations  can  now  be  made.  First  of  all,  it  is  easy 

to  see  from  (Eq.  30)  that  the  limits  obtained  in  this  fashion  are  stricter  • 

than  the  ones  for  convergence  of  the  mean  of  k.  .  Furthermore,  while 

l  ,n 

previous  results  showed  a  dependence  of  the  upper  limit  on  E^_^  only,  and 
thus  the  upper  limit  for  a  could  be  made  independent  of  the  position  in 
the  lattice  by  normalizing  it  with  respect  to  E^_^  [31,  Equation  (30) 
shows  that  a  is  also  dependent  on  the  crosscovariance  of  f^_^(»)  an<* 

bi-l<*>* 

Thus,  in  order  to  make  the  limits  for  the  stepsize  independent  on 
the  position  in  the  ladder,  normalizing  with  respect  to  the  right-hand 
side  of  (Eq.  30)  requires  the  extra  computation  of  C^_^.  As  is 

bounded  by  the  Schwarz  inequality  as: 

« 

0  -  ici-il  s  Ei-i  <31) 

l 

two  extreme  cases  can  be  considered.  If  is  very  small  (input 
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sequence  extremely  random)  the  upper  bound  for  a  in  (Eq.  30)  becomes 
approximately  equal  to  •=— —  ,  while  if  the  input  sequence  is  almost 

El-1  2  ! 

deterministic  (singular  process),  the  limit  for  a  approaches  •=■  — - . 

3  Ei-1 

This  indicates  that  faster  convergence  can  be  obtained  for  lower  signal- 
to-noise  ratios,  a  result  which  seems  to  be  supported  by  some  experimental 
work. 

3.  The  Misadjustment 

The  misadjustment  due  to  the  ith  stage  in  the  lattice  filter  is 
defined  by  (Eq.  21).  Since  (Eq.  23)  shows  that 

E{lki.»-kil2)Ei-i  (32> 


it  is  clear  that  the  normalized  misadjustment  is  nothing  but  the  norm  in 
the  fixed-point  problem  (Eq.  26): 


2  A 


M‘ 


m. 


"i-1 


lim 

n-*» 


(33) 


The  remark  can  be  made  that  the  fixed-point  of  T^,  the  "ideal"  contra¬ 
ction  mapping  given  by 

T£  k  ■  [1-  a  Ei_1]k  -  a[Ci_1  +  k±  Ei_1]  (34) 

has  as  (only)  fixed-point  the  zero  element  because 

Tt  0  -  0-a[Ci_1  +  ki  Er] 

-  0  (35) 

by  the  definition  of  k^  (Eq.  16).  Thus,  zero  misadjustment  can  be  obtained 

A  A 

if  the  estimates  in  and  converge  to  the  actual  values. 


i 


16- 


December  8,  1982 
GRLS :LHS :cac 


This  will  be  Che  case  if  a  different  estimate  Is  used  as: 


(36) 


and 


ci_l(n) 


n  A  Ifi-l(j)bi-l(J)] 


(37) 


and  if  the  process  e(*)  is  correlation  ergodic  [18].  Note  that  in  this 
case,  ^  satisfies  the  established  conditions  for  stochastic  approximation 
[19]. 

In  order  to  compute  the  misadjustment  of  the  algorithm  (Eq.  17), 
the  norm  of  (Eq.  26)  is  needed. 


lki,n+l’kiH  l^i,n  <ki,n“ki*  H 


(38) 


This  norm  is  computed  in  Appendix  B  and  the  (normalized)  misadjustment 
is  found  to  be: 


2 

‘i 


2  2 
Ei-1 ~  Ci-1 

Ei-1 


-i2 


(39) 


This  expression  indicates  that,  while  the  limits  on  a  can  be  made  in¬ 
dependent  of  the  stage  i  by  appropriate  normalization  of  the  stepsize, 
such  normalization  will  not  result  in  a  stage-independent  misadjustment. 


Again,  two  extreme  cases  occur.  If  the  input  sequence  is  extremely  random 

a  E.  , 

( I  ci_i  I  <<Ei_i^ »  *Eq'  wil1  equal  approximately  j  ,  while  if 

the  input  sequence  is  almost  deterministic,  the  misadjustment  approaches 


zero.  Therefore,  is  bounded  by: 
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2  a  ®i-l 

0  <  Mi  <  2— cTe^' 


In  this  result,  ^  is  the  average  power  in  steady  state  conditions 
at  the  input  of  the  ith  stage.  For  the  usual  case  where  a  is  normalized 
with  respect  to  the  upper  bound  in  (Eq.  40)  becomes 

0  ‘  Mi  <  20^0  Ei-1 

In  order  to  compute  the  total  misadjustment  at  the  output  of  the 

cascade,  (Eq.  13)  can  be  used.  Essentially,  (Eq.  13)  corresponds  to  the 

model  shown  in  Figure  3.  The  total  output  power  in  Figure  3  is  bounded 
2 

by  (Eq.  13)  where  and  k^  are  (since  the  lattice  converged)  equal  to 
2 

m^  and  k^,  respectively: 

»  t«J+(l-|k1l)2)E0<  Ej,  <  [mJ+(l+|ki|)2jEo 

Since  the  optimal  output  power  is  given  in  (Eq.  8b),  the  output  mis¬ 
adjustment  normalized  with  respect  to  the  total  output  power  is  given  by: 


N  F m2+  (1-  | kt I ) 2 


N  [V+U+lkJ)2 


If  the  input  sequence  is  very  random,  |k^|<<  1  and  both  bounds  approach 
the  same  value.  In  this  case,  we  will  have  approximately 

„  «  2 

ra  *  it  (1  +  m.  )  -  1 

i-1  1 

If  a  is  normalized  with  respect  to  E^,  which  is  the  case  in  most 
applications,  we  have  from  (Eq.  41): 


i|MM| 1 

'!  ••  j 
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2  ct 
mi  “  2(1  -  a) 

and  finally 


Ungerboeck  [21]  and  Widrow  [22]  showed  that  for  the  LMS  algorithm, 
the  misadjustment  is  linearly  dependent  on  the  order  N.  The  exponential 
dependence  on  N  for  the  lattice  suggests  that.  In  general,  a  higher  mis- 
adjustment  can  be  expected  for  the  lattice  filter. 


4.  The  Convergence  Rate 

In  order  to  compute  the  convergence  rate  of  the  ith  adaptive  element 
the  homogeneous  part  of  the  state  equation  in  (Appendix  B  -  Eq.  B8)  has 
to  be  considered.  This  is  denoted  by: 

fA  -  2  o  C  1 


x(n  +  l)  * 


S(n) 


1  -  a  E 


i-1 


5  M  x(n) 

The  convergence  rate  is  determined  by  the  smallest  eigenvalue  of  the 
matrix  M  [20].  Since  it  is  straightforward  that  this  smallest  eigen¬ 
value  is  equal  to  A,  the  "slowest"  mode,  which  determines  convergence  of 
the  squared-error,  is  given  by: 

« n 

or  *  A 
i,n 

Substituting  for  A  using  (Appendix  A  -  Eq.  A3)  and  assuming  small  a, 
an  approximate  expression  for  m  is  obtained: 
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Comparing  (Eq.  49)  Co  the  expression: 


Che  time-constant  of  the  adaptation  is  given  by: 


2cE1_1-c.2[2eJ_1*C2_1] 


(49) 


(50) 


In  the  extreme  case  of  almost  deterministic  inputs  (high  signal-to-noise 
ratio),  this  becomes: 

= - H -  (51) 

2aEi_1[l--faE1 

and  in  case  of  highly  uncorrelated  input  sequences  (low  S/N  ratio), 
becomes: 

T  , - 2 -  (52) 

2aEj.jIl-aEj.jI 

Thus,  a  slower  adaptation  is  obtained  for  deterministic  sequences  than 
for  uncorrelated  sequences. 
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CONCLUSION 

Our  results,  which  are  based  on  a  mean-square  convergence  criterion, 
suggest  that  the  usual  approach,  which  consists  of  the  study  of  conver¬ 
gence  properties  of  the  mean  of  the  reflection  coefficients,  may  not  be 
sufficient  for  a  number  of  practical  applications.  The  mean-square 
approach  shows  that  convergence  is  only  obtained  under  stricter  conditions 
than  previous  studies  suggested.  The  convergence  rate  of  the  adaptive 
gradient  lattice  is  better  for  highly  uncorrelated  input  signals  (low 
signal-to-noise  ratio)  than  for  almost  deterministic  sequences  (high 
signal-to-nolse  ratio),  but  a  higher  misadjustment  can  be  expected  if  the 
signal-to-noise  ratio  is  low.  The  total  output  misadjustment  varies 
exponentially  with  the  filter  order.  Since  the  output  misadjustment  for 
the  basic  LMS  algorithm  varies  linearly  with  the  number  of  taps,  the  excess 
mean-square  error  for  the  lattice  can  be  expected  to  be  larger,  especially 
for  higher  order  filters. 


> 


% 
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APPENDIX  A 


Note  that: 


E{  1 1  -  ct  Ei_1(n)  |  2} 

=  E{l+a2  j Ei_1  <n)  | 2  -  2  aE^^n)} 

=  1+a2  E^E^^n)  j2}  -  2aE1_1  (Al) 

Now,  using  the  well-known  expression  for  4-th  order  Gaussian  moments  [17], 
this  becomes  from  (Eq.  18) : 


E{|ei.1(n)|2}  -  i  Et  I  ^(11)+ 8^(0)  I2} 

=  i[E{fi-l(n)+bi-l(n)} 

+  2  E{f2_1(n)  b2_1(n)}J 

■i[3ELi+3E?-i+2Ei-i+4Ci-i] 


(A2) 


Substituting  this  in  (Appendix  A  -  Eq.  Al),  we  finally  obtain: 
E{  1 1  -  a  E.^Cn)  | 2}  -  a2  [2  E2^  +  C2^]  +  1  -  2  a  E^ 


(A3) 
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Using  (Eq.  25),  the  following  recursion  can  be  written: 

II  ki,„+i  -  kJI 2  ■ E  <  1 11  -  “  8i-i (n)  i  (ki -  ki> <»> + ki  h-i 

‘  EfJtl-aE^U)]2  (ki  n-ki)|2} 

+  a2E{[ew(n)+k1lw(n)]2} 

-  2aE(tl-aEi_1(n)]  [C^n)  +  k±  ^ ^(n)  ]  -k^} 


(n)]|2} 


(Bl) 


If  we  assume  that  the  repetition  rate  of  the  algorithm  is  small  enough 

such  that  (k.  -k.)  is  uncorrelated  with  C.  , (n)  and  E.  , (n) ,  (Appendix  B  - 

i,n  l  l-l  l-l 

Eq.  Bl)  will  become: 

II ' kl,„+i  -  kt H : 2  * E<  I  ■ 1  -  “ : I ' 2 i  II ■ ki,„  -  ki II : 2 

+  at2  E  {  [2i.1(n)  +  k£  §i_1(n)  ]2  } 

-  2  a  E{  [1  -  a  ^(n)  ]  [C^n)  +  k±  E^n)  ]  } 


E{kn  n'M 

i,n  l 


(B2) 


The  coefficients  in  (Appendix  B  -  Eq.  B2)  can  again  be  computed 
using  the  expression  for  4-th  order  moments  of  Gaussian  processes: 


EUVi(n)+kiE1_1(n)]2} 

-  E(|ei_1(n)|2}+k2  E{|Ei  1(n)|2}  +  2Vi£{Ci ^(n)  E^n)} 


(B3) 


In  (Appendix  B  -  Eq.  B3)  we  have: 


7 
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E(  |Ci_1(n)  |2}  -  Etf^n)  b^n)} 


2  2 
Ei-1  +  2  Ci-1 


E^n)  Ei_1(n)}  =|-E{fi_1(n)  bi_1(n)  f^Cn)} 


+  -|E{fi_i(n)  bi_1(n)  bi_1(n)} 


3C1-1  Ei-1 


Using  (Eq.  B4)  and  (Eq.  B5)  and  (Eq.  A2) ,  (Eq.  B3)  becomes: 


EUC^OO+lCi  Ei_1(n)]^} 


'  Eli+2Ci-i+2ki  Ei-i+kJ  cli 


+  6  k±  Ci_1  E1_1 


and  with  the  value  for  the  optimal  reflection  coefficient  from  (Eq.  16), 
(Appendix  B  -  Eq.  B6)  reduces  to: 


E{tCi_1(n)+kiE1_1(n)}2  }-  E2 ^  +  2  C2^  +  2  C2  x  +  -|=^  -  6  C2_x 

Ei-1 

C4 

=  e2  _2c2  +  fi=l 

i-1  i  i-1  2 

Ei-1 

rc2  _r2  “I2 


“  r2  -  r2 

Ei-i  ci-i 

Ei-i 
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The  recursion  in  (Eq.  B2)  can  therefore  be  written  as: 

llki,n+l'kiH2“  A  II ki ,n  ~ ki II 2  +  01,2  ® 

-  2  a  C  E{k  -  k. }  (B8) 

i,n  i 

where  A  is  given  by  (Eq.  A2) ,  and  B  is  given  by  (Eq.  B7) .  Equation  (B8) 
can  now  be  considered  as  a  state  equation  with  state  variable  given  by: 

(B9) 


x(n)  « 


E<kl,„-ki) 


Transfer  analysis  of  this  state  equation  [20]  and  application  of  the 
final-value  theorem  of  the  Z-transform  lead  to  an  expression  for  the 
(normalized)  misadjustment  due  to  the  ith  stage: 


2  a2B 

mi  =  I^A 


(BIO) 
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